Solitons and thermal fluctuations in strongly non-linear solids 
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We study a chain of anharmonic springs with tunable power law interactions as a minimal model to 
explore the propagation of strongly non-linear solitary wave excitations in a background of thermal 
fluctuations. By treating the solitary waves as quasi-particles, we derive an effective Langevin 
equation and obtain their damping rate and thermal diffusion. These analytical findings compare 
favorably against numerical results from a Langevin dynamic simulation. In our chains composed 
of two sided non-linear springs, we report the existence of an expansion solitary wave (anti-soliton) 
in addition to the compressive solitary waves observed for non-cohesive macroscopic particles. 
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In linear elastic solids, phonons are the basic mechan- 
ical excitations responsible for energy propagation. By 
contrast, an aggregate of macroscopic grains just in con- 
tact with their nearest neighbours constitute a novel 
elastic material where solitary waves or shocks replace 
phonons as the basic excitations [ll-Q- The origin of these 
strongly non-linear waves can be traced to the fact that, 
unlike the case of harmonic springs, the repulsive force 
between two grains in contact does not depend linearly 
on the relative compression. So far, little effort has been 
directed to determine the fate of these strongly non-linear 
excitations in a background of thermal fluctuations be- 
cause temperature is clearly not a parameter relevant to 
the elastic response of macroscopic grains. 

However, granular aggregates at zero pressure are just 
one example of a broader class of materials that can be 
prepared in a unique mechanical state called sonic vac- 
uum [l|. This term originally coined by Nesterenko in 
the context of strongly non-linear granular chains des- 
ignate a material characterized by a vanishing elastic 
moduli and linear speed of sound [l|, Ij-Q . Grafted col- 
loidal particles Q and ultra-cold atoms in optical lattices 
[10| are microscopic systems that allow for tunable non- 
linear interactions, while being naturally coupled to a 
source of fluctuation (thermal or quantum) . These fluc- 
tuations restore rigidity and generate long wavelength 
phonon modes [li|,[l2l- However, the physics of very high 
amplitude strain propagation is still predominantly non- 
linear and resembles the state of sonic vacuum perturbed 
by background fluctuations. This extreme regime is par- 
ticularly relevant for energy transport in some biological 
systems, where transport occurs through localized non- 
linear excitations with energy significantly higher than 
the thermal energy [1^ [ij] . 

Moreover, systems such as polymer networks and col- 
loidal glasses undergoing an unjamming transition arc 
also characterized by vanishing elastic moduli as the coor- 
dination number or packing fraction are lowered towards 
the critical point [l5| . The effect of thermal fluctuations 
on the non-linear response of materials undergoing an 
unjamming transition is relatively unexplored, despite 
they are obvious examples of a sonic vacuum state at 



zero temperature |16l - ll8| . Note, that in the case of jam- 
ming the linear elastic moduli can be lowered towards 
zero even if the microscopic interactions arc harmonic, 
simply because there arc not enough forces to prevent 
floppy motions. 

In this article, we focus on strongly non-linear mechan- 
ical waves propagating in a background of small thermal 
fluctuations, a non-equilibrium problem that lies outside 
the realm of perturbation theory. The starting point 
of conventional perturbation methods is a linear elastic 
solid, possibly at finite temperature, perturbed by small 
anharmonic terms. By contrast, we adopt as a starting 
point the fully non-linear state of sonic vacuum whose 
elementary excitations arc long-lived solitary waves [19| . 
Subsequently we switch on temperature as a small per- 
turbation that creates a background of thermal fluctua- 
tions. 

As a minimal model that is analytically tractable, we 
study impulse propagation in a one dimensional lattice 
of non-linear springs with a tune-able power law interac- 
tion. In addition to the compressive solitary waves seen 
in a lattice of macroscopic grains with one-sided repul- 
sive interaction, we report an accompanying anti-solitary 
wave solution for the lattice of non-linear springs with 
two sided interactions. By coupling the lattice to a heat 
bath, we then study the effects of the thermal fluctua- 
tions on the leading solitary wave generated in response 
to an impulse of energy much higher than the background 
thermal energy. Our approach in a nutshell is to treat the 
solitary wave as a quasi-particlc and derive an effective 
Langevin equation that describes its stochastic dynam- 
ics. We corroborate our analytical predictions for the 
damping rate and thermal diffusion of the solitary waves 
with Langevin dynamic simulations. 

Model Consider flrst, a lattice of beads of mass to, 
where nearest neighbours interact with a one sided non- 
linear potential of the form V{S) = —S", with S being 
the overlap between neighbouring beads, k the force con- 
stant and a > 2 the exponent of the non-linear potential 
(a very large may be used to obtain the limit of hard 
sphere interaction). If the initial overlap; (5 — > 0, the ef- 
fective spring constant defined as the second derivative of 
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FIG. 1: left:{a.) Velocity profile of the compressive solitary 
wave (SW) generated in an athermal chain of beads with 
one sided interaction, (b-c) Velocity profiles showing the for- 
mation of a train of SW-ASW pair for two-sided non-linear 
springs. A single particle is initially given an impulse to the 
right, generating a train led by a SW moving in the direc- 
tion of the impulse (b), while simulatenously generating a 
symmetric train led by an ASW moving in the opposite di- 
rection, right: The energy momentum relation for the leading 
SW/ASW in the three cases shown in the left panel following 
the energy (E) momentum (P) relation E — 
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inset: 



Zoom in of the leading SW-ASW pair from left panel (b). 



the potential vanishes, and as result, the lattice does not 
support linear sound at zero temperature. Consequently, 
mechanical strains generated in response to an impulse 
evolve into long-lived non-linear solitary waves and the 
solitary wave energy E and momentum P [2^ satisfy the 
classical relation E 



where rricff is the mass of the 



solitary wave quasiparticle [1|, |7| . 

In Fig. ([T]), we demonstrate that the compressional soli- 
tary wave (SW) excitation discovered by Nesterenko in 
a chain of non-cohesive beads is also seen in a lattice 
of springs with two sided interactions. However, unlike 
the case of a one sided potential, each compressive SW 
generated in response to an impulse is accompanied by a 
corresponding expansion solitary wave (formed by local 
stretching of springs) of the same magnitude but moving 
in the opposite direction. This anti solitary wave (ASW) 
is not sustained by beads interacting with purely repul- 
sive potentials - the beads would merely loose contact. 

In Fig. dl]) left panel, we show the SW/ASW excita- 
tions for (a)beads , (b-c) particles connected by springs. 
In Fig. ([T]), right panel, we plot the energy- momentum 
relation for the leading SW/ASW demonstrating that 
SW excitations in a lattice of repulsive beads have the 
same effective mass nicff as a SW and ASW in two-sided 
springs. As shown in Fig. ([1]) left panel (b-c), the leading 
SW-ASW generated in response to an impulse imparted 
to one of the particles towards the right (direction of ar- 
row) is followed by a train of alternating SW-ASW's exci- 
tations, of progressively smaller magnitudes. The smaller 
SW/ASW's are generated as the particle initially im- 



parted the impulse, recoils with its left-over energy. This 
process is repeated several times, leading to the genera- 
tion of the train of smaller excitations. Since the speed 
of propagation depends upon the amplitude, the SW and 
ASW that start propagating together initially (appearing 
bounded), eventually separate and become clearly distin- 
guishable. 

Langevin Equation The classical energy-momentum re- 
lation (see Fig. ([T|) right panel) satisfied by the SW mo- 
tivates the interpretation of the solitary wave as a quasi- 
particle [H, Q- For small perturbations, the SW can still 
be treated as a quasi-particle provided the effects of the 
perturbations accrue gradually such that the SW retains 
its functional form. We now apply this adiabatic approx- 
imation to derive an effective Langevin equation for the 
SW quasi-particle when the lattice of springs is coupled 
to a heat bath. Recall first, the Langevin equation for a 
particle of mass m undergoing Brownian motion in one 
dimension is 
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Here, £■, K are the total and kinetic energies respectively, 
C, /3 are the dissipation and diffusion coefficients related 
via the fluctuation dissipation theorem /3^ = 2(kBT/m^, 
where ks is the Boltzmann constant. N{0, 1) is a normal 
random variable with mean and variance 1 , and encap- 
sulates the effects of random fluctuations during the time 
interval t,t + dt. For a free particle of unit mass moving 
with speed v, E = K = ^v^ and upon substituting in 
Eq. ([T|), we recover the Langevin's equation convention- 
ally expressed as the rate of change of momentum of the 
particle j23 |. 

We now derive an equation analogous to Eq. ( [1]) for 
the compressive solitary wave quasi-particle. Identifying 
the lattice spacing a as a characteristic length scale and 

Lo = a/;^^""^ as an inverse time scale, the equation of 

motion for the compressive displacement field 4'{x, t) in 
dimensionless units reads- 
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where subscripts denote partial derivatives with respect 
to space X and time t. Eq. ([3]) is a simplified form of the 
Nesterenko equation [l|, y] . The first two terms express 
the rate of change of momentum while the third term 
represents the force. Although the solitary wave solution 
to Eq. ^ is not exact (lacking compact support), Eq. 
([2]) provides a good approximation while being analyt- 
ically more tractable especially since we are interested 
in keeping the non- linear exponent a general [2, l3| . Note 
that the equation for the ASW (stretching) is obtained by 
modifying the third term +[(-0x)"~ ]x -^ ~[{(i>xT~ ]x 
in Eq. ^. Upon substituting 5 — —(px for the compres- 



sive SW OT 6 = (px for the expansive ASW, wc find the 
same functional forms for the solitary wave solutions in 
both cases. 

In analogy with the Langcvin equation for a particle, 
we model the coupling to a heat bath as the sum of two 
contributions- an external drag and a random fluctuating 
force, phenonienologically introduced into the equation of 
motion as : 

<j)tt - J^(j)xxtt + [{-(l)xT~^]x = -7 ( f/** - Y^0tex j + 



-J^ L{x, t; t + dt)- -^Vx{x, t- t + dt)] (3) 



where 7 = -^ is the dimensionless drag coefficient that 
couples to the momentum 11 = (0t — j24'txx)- It is use- 
as the ratio 



ful to define a coupling constant F — ^^ j, 
of potential to thermal energy in terms of which, the 
dimensionless diffusion coefficient is I? = ^. The last 
(noise) term on the right of Eq. ^ in conjunction with H, 
satisfies the fluctuation dissipation theorem [23|. Here, 



r](x, t]t + dt) is a Gaussian random noise during the time 
interval t,t + dt with the moments (77(2;, t;t + dt)) = 
and {r]ix,t;t + dt)r]{x',t';t' + dt')) = S{x~ x')5{t-t') re- 
spectively, where angular brackets denotes ensemble av- 
eraging. 

To study the propagation of the SW in a background of 
thermal fluctuations, we now make a working assumption 
based on the quasi-particle approximation to the SW: 
whenever the energy of the SW, E = -Bsw ^ ksT, the 
SW functional form is unaltered and only its amplitude 
A{t) becomes time-dependent. The amplitude A{t) is 
the collective variable for the SW quasiparticle and other 
properties of the solitary wave, such as its energy and 
momentum may be determined from it. Note, the width 
of the SW is independent of its amplitude and therefore 
we do not consider its time dependence |23j . 

From Eq. ([3]), the conserved energy is 



E ^ I dx -> 



-L</)L + -(-0.)", (4) 

24 a 



and the energy of the SW may be obtained by integrating 
Eq. (|4]) over the width of the SW of order W. (This avoids 
including the energy of small SW that separate from the 
main wave) . Using Eq. Q , the rate of change of energy 
is (see SI for details). 
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where, K is the kinetic part of the energy. 
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The last two terms on the right of Eq. ([5]) describe the 



possible mechanisms of decay of the SW by "friction" 
from the heat bath (7) and the "phonon drag" induced by 
the thermal motion of the chain a{t). The first term is the 
fluctuating part of the energy. In the following, we make 
the assumption (verified numerically) that the coupling 
to the heat bath is more important and therefore, ignore 
a{t) (see SI for details). 

Solving for the SW solution from Eq. ^ , we find (pt ~ 
Alp, where ip{x,t) = sech°^^ (^TF^) i^ ^^^'^ functional 
form of the SW with amplitude A, speed Vs = A~^ 



and width W 



in units of the lattice spacing. 



\/3(a-2) 

The SW energy, kinetic energy, and momentum may now 
be expressed in terms of the collective variable A. The 
energy \^ E — J^ ^ dx 4''^{x, t) = A^Ie, and from the the 
virial theorem, K = -^E = -^IeA'^. Additionally, 

' a-t-2 a-\-2 ^ -^ ' 

the momentum P ~ A J dx ip{x,t) = Alp, where Ie-Jp 
are constants obtained by integrating over all space [3 . 
Substituting for E and K in terms of A allow us to cast 
Eq. ([5]) into the form of an ordinary Langevin equation 
(with additive noise) for the collective variable A{t), 
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where, 4> = (pi^^t)- EQ- © is the central result of our 
work whose analytical predictions we derive and test nu- 
merically in the next sections. The first term can be 
written as -y=-qA{t,t + dt), where r}A is a white noise sig- 
nal; that is, its correlations are given by {r]A{t)rjA{t)) = 
(a+2')i r ^^^ ~ *)■ Using the fact that the correlations of 
•q{x,t) are described by delta-functions, the correlations 
of rjA{t) can be related to the kinetic energy Eq. ([6|), 
which can be replaced by -^^^lEAit)'^ . 

Mean and variance Taking the expectation value of Eq. 
([7]), we find the solitary wave amplitude decays as 



{A) = Aoe- 



(8) 



where, Ao is the initial solitary wave amplitude. Here, 
owing to the noise term r]{x, t; t+dt) (which acts between 
times t;t + dt) and (j)t{x,t) (which is a solution at time 
t) being statistically independent, the expectation value 
{il{x,t;t + dt)<pt{x,t)) = 0. Note, the effective damping 
rate is independent of inverse temperature F but rescales 
with the exponent of the non-linear potential a. 

Similarly, we solve for the variance of the solitary wave 
amplitude or equivalently, the variance in the square root 
of energy to find (see SI for details) 



var(vi?) = 



2ar V 
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(9) 



The coefficient becomes !^^ when the energy is not mea- 
sured in units of fca", so this expression is analogous to 
the velocity variance of a Brownian particle. Note, for 
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FIG. 2: left: Snapshot of two solitary waves in a background 
of thermal fluctuations for a = 2.5 (red) and a — 2.2 (black). 
right: the attenuation of the solitary wave as a function of 
time for various values of 7,r, and a. When not indicated, 
a — 2.5, r = 10* and 7 = 5 x 10~5. The data for these values 
is shown with green circles; the other data corresponds to 
changing one parameter, a = 2.2 (black circles), F — 1.6x10'', 
(brown squares) and 7 = 7x10"'' (red circles) compared with 
the analytic expression in Eq. (|8]) represented by solid lines. 
The initial Esw = 0.5. 



FIG. 3: left: the linear increase in the variance of the solitary- 
wave square root of energy, as a function of time for various 
values of a, 7, and F. Where not indicated, a = 2.5, F — lO'' 
and 7 — 10^*. right: The collapse of the numerical data 
spanning F — {.06,0.1, lIxlO** for the variance of the square 
root of the solitary wave energy plotted as function of time on 
a linear scale and compared against the analytic expression 
in Eq. (fTU)) . Where not indicated, a = 2.5, F = lO" and 7 = 
5x10^''. The initial Esw = 0.5. For each curve, a different 
value of the background noise has to be subtracted. 



large a, the SW is effectively one particle wide and thus 
Eq. ^ captures the correct thermal equilibration of the 
particle energy with the heat bath. However, for the dy- 
namics of the SW, Eq. ((9]) is only useful as long as the 
SW is identifiable against the background thermal en- 
ergy, that is -Esw ^ T~^. We thus linearize Eq. © to 
find how the variance accumulates at early times and use 
this expression to compare against numerical data, 
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(10) 



Simulations We consider a one dimensional chain con- 
sisting of A^ ~ 1024 particles each having a mass m placed 
regularly on a lattice with spacing a (spring rest length) 
interacting pair-wise with a nearest neighbour interaction 
V{6) = — ((5)", where S is the compression/stretching in- 
duced during the dynamics. We model the coupling to 
the heat bath by numerically integrating Eq. ( [T]) for each 
particle using the predictor-corrector algorithm |2y|. In 
thermal equilibrium the mean kinetic energy is KE^ ■-^— 

and potential energy is PE~ ^^ , where their ratio sat- 
isfies the virial relation, see Supplementary Information 
(SI), Fig. (1), inset for a = 5/2. In the following, all 
numerical data is presented in dimensionless units, en- 
semble averaged over 1000 samples. 

In SI, Fig. (1), we plot the sound speed from the slope 
of the dispersion curves for a = 5/2 for a range of F. At 
thermal equilibrium, the mean kinetic energy and hence 
the temperature T satisfy the virial relation T 
where St is the average displacement of the particles in 
duced by thermal fluctuations. Defining the sound speed 



Urp, 



c as the second derivative of the induced potential energy 
leads to the relation, c^ ^ T~^ [12]. For a = |, we find 

c ^ Ft^ ^ {ksTy^^'^, closely matching the linear fit in SI 
Fig. ( 1). Thus, coupling the lattice of non-linear springs 
that is initially in its state of sonic vacuum (implying 
the absence of linear sound) to a heat bath, leads to hy- 
drodynamical sound modes with a linear sound speed 
that scales with the temperature of the heat bath |13 |. 
Note, setting a = 2 (harmonic springs) yields a sound 
speed that is independent of temperature while the limit 
a — >■ cx) yields c ^ (fc^T)^'^, a result in agreement with 
the entropic elasticity for hard sphere colloidal crystals 

El. 

Once the lattice reaches thermal equilibrium, we excite 
a solitary wave (SW) by imparting one of the particles an 
initial energy of order i?sw = 0.5 in dimensionless units. 
In Fig. ([2]) left panel, we show a snapshot of two SWs at 
the same time, propagating in a background of thermal 
fluctuations for a = 2.5 (red) and a = 2.2 (black). We 
see that the SW with lower a is wider and moves faster 
for the given amplitude, in agreement with the analytic 



widths W 



and speeds Vs ^ A~ 



In Fig. 



%/3(a-2) 

([2]) right panel, we plot the numerical data (symbols) for 
the attenuation of the SW amplitude as a function of 
time for various values of 7, F and a and we find a very 
good match to the analytic expression in Eq. ([5]) (solid 
lines). For the range of F explored, we find the damping 
rate is independent of temperature F but depends on the 
environmental drag 7 and a. 

In Fig. ([3]) left panel, we show the linear increase in 
the variance of SW amplitude (or the square root of its 
energy) as a function of time for multiple values of a , 



7 and F. The right panel compares the numerical data 
to Ea. ([TU)) . After rescaling by the theoretical value of 
the slope and shifting each curve, the data collapses on a 
line. These predictions apply when 7^ ^ 1 and provided 
that the SW energy is much greater than the thermal 
energy T~^ [231 ■ Furthermore, the mean damping rate 
and variance are independent of the initial SW energy. 

To conclude, we find that a lattice of two sided non- 
linear springs generates a pair of solitary and anti-solitary 
waves in response to an impulse. By coupling the lattice 
to a heat bath, we study the propagation of a solitary 
wave with energy much greater than the thermal energy 
by deriving an effective Langevin equation for the solitary 
waves. We calculated the damping rate and the growth 
rate of energy fluctuations, and verified our results with 
simulations. 



SUPPLEMENTARY INFORMATION 



Derivation of energy equation and solitary wave ampli- 
tude variance We first derive the variance of the ampli- 
tude of a soliton from Eq. (5) in the main text, and then 
provide a rationale for the phonon induced drag <y{t). 
Re-defining, D = ^^^p , we solve for the variance in the 
solitary wave amplitude by first evaluating the differen- 
tial ^[yl^] ^ A^{t+dt)-A^{t), by substituting A(t-hdt) == 

^ (1 - £+2'^*) +VDdiJdx r){x, t- 1 + dt) (0t + -^<Pxt) 
from Eq. (7) in the main text and retaining terms to 
order 0(dt)[23,[2l, 
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dxdx' ri{x,t;t + dt)ri{x' ,t;t + dt) ((j)t{x,t) + —-=(f>tx{x,t) j i 4>t{x' ,t) + —-=(j)tx'{x',t) 
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Taking the expectation value, the second term on the 
right vanishes (as discussed for the mean) and using the 
property that the noise term is delta-correlated in space, 
we obtain 



Ddt / dx 



1 



d[{A')] = -^{A')dt 



The last term when expanded gives twice the soli- 
tary wave kinetic energy 2K(Eq. (6)) plus an integral 
J dx(j)t4>tx, that vanishes by symmetry for the SW 
Moreover, the SW kinetic energy is related to 
its total energy via the virial relation K = -^E. Hence, 
we obtain the ordinary differential equation correct to or- 
der dt - 
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Solving, the differential equation subject to the initial 
condition {A'^)t=o = Aq and substituting for D, we ob- 
tain. 



(A^) = Ale'^' + 
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2/_EaF 
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(12) 



The energy of the soliton changes because of the inter- 
actions with the noise, friction, and conduction from the 
rest of the particles in the chain. The energy is 
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x(t)+l/2 



x(t)-l/2 
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^</)L + ^(-0.)"dx, (13) 



where / is the width of the soliton (which can be arbi- 
trary as long as it is wide enough to include the whole 
soliton-making I larger just adds a small offset to the 
energy of the soliton by including the thermal fluctua- 
tions of the particles in the background). The deriva- 
tive with respect to time has one contribution from dif- 
ferentiating inside the integral. If one evaluates this 
derivative and integrates by parts, the integrand becomes 
'/'(('/'tt - T^4>uxx - {-fpxT'^)- Next substitute the differ- 
ential equation for (j) into Eq. (3) above and integrate by 
parts to get the first two terms from Eq. (5) in the main 
text. The boundary terms will be part of (j{t). 

Due to the translation of the soliton, we pick up an- 
other contribution to ^. After the interval moves it 
extends more on the right and less on the left (the end- 
points move the distance Vgdt if 1;^ is the speed of the 
soliton). This changes the energy by Vsdt{^<j)l + j^(j>^^ + 



terms gives <T(t), 



Combining this with the boundary 
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Fluctuation induced rigidity To extract the equilib- 
rium properties in the thermalized state, we define the 
longitudinal current density of particles as j{x,t) = 
Vi(t)S{x — Xi{t)), and its Fourier transform 
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where k is the longitudi- 
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FIG. 4: The sound speed computed from the dispersion curves 
for a range of F (inverse temperature) for q = | . The circles 
are from numerical simulation while the dashed blue line is a 
linear fit, giving a slope of 0.11 for a — ^, close to the expected 
value ^^. The inset shows the kinetic and potential energy 
approaching thermal equilibrium, where their ratio satisfies 
the virial relation KE/PE = a/2. 



nal collective mode along the a;-direction. Thus, the cor- 
responding longitudinal current density auto-correlation 
functions is C{k,t) = {j*{k,0)j{k,t)), where the an- 
gular brackets denote ensemble averaging over the ini- 
tial time. The longitudinal power spectral density is 
then obtained as the Fourier transform of the respective 
current density auto-correlation functions as, P{k,Lu) = 
/_ dt e'^'^*C{k,t). The Fourier transforms defined above 
are evaluated using fast Fourier transform from simula- 
tion data. The sound speeds in SI Fig. (j4|) correspond 
to the linear part of the dispersion curves, obtained by 
projecting the power spectral densities on the frequency 
(a;)- wavenumber (k) plane. 



Acknowledgments NU acknowledges financial sup- 
port from FOM. 



[1] V. F. Nesterenko, J. Appl. Mech. Tech. Phys. 5, 733 

(1984). 
[2] P. Rosenau, Phys. Lett. A 118, 5 (1986). 
[3] L. R. Gomez, A. M. Turner, M. van Hecke and V. Vitelh, 

Phys. Rev. Lett. 108, 058001 (2012). 
[4] V. F. Nesterenko, Dynamics of Heterogeneous Materials 

( Springer- Verlag, New York, 2001). 
[5] C. Coste, E. Falcon, and S. Fauve, Phys. Rev.E 56, 6104 

(1997). 
[6] C. Daraio, V. F. Nesterenko, E. B. Herbold, and S. Jin, 

Phys. Rev. E 72, 016603 (2005). 
[7] S. Job, F. Melo, A. Sokolow, and S. Sen, Granular Matter 

10, 13 (2007). 
[8] S. Sen, J. Hong, J. Bang, E. Avalos, and R. Doney, 

Physics Reports 462, 21 (2008). 
[9] A. Yethiraj, Royal Soceity of Chemistry 3, 1099-1115 

(2007). 
[10] D.E. Chang, J.I. Cirac and H.J. Kimble, Phys. Rev. Lett. 

110, 113606 (2013). 
[11] Z. Cheng, J. Shu, W.B. Russel and P.M. Chaikin, Phys. 

Rev. Lett. 85, 7 (2000). 
[12] O.V. Zhirov, A.S. Pikovsky and D.L. Shepelyansky, Phys. 

Rev. E 83, 016202 (2011). 
[13] A.S. Davydov, J. Theor. Bio. 66, 379 (1977). 
[14] T. Dauxois and M. Peyrard, Physics of solitons (Cam- 
bridge University Press, 2006). 



[is; 

[16 
[17 
[18 
[19 

[20 

[21 
[22; 
[23; 
[24 
[25; 
[26; 
[27 



Chem. Phys. 138 
R. Nagel, Phys. 



C. S. O'Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, 
Phys. Rev. E 68 011306 (2003). 

M. Sheinman, C.P. Broedersz and F.C. MacKintosh, 
Phys. Rev. Lett. 109, 238101 (2012). 
I. Atsushi, L. Berthier and G. Biroli J. 
12A507 (2013) 

N. Xu, T. K. Haxton, A. J. Liu, and S 
Rev. Lett. 103 245701 (2009). 

We are using "soliton" and "solitary wave" interchange- 
ably to refer to a localized lump of of energy; we do not 
assume anything about the equations being integrable. 
For many waves or solitons, differentiating E with respect 
to P would give the group velocity or the speed of the 
wave's propagation, but this does not work for this case. 
{P is not the canonical momentum.) 
M. Manjunath, A. P. Awasthi, and P.H. Geubelle, Phys. 
Rev.E 85, 031308 (2012). 

N.G. Van Kampen, Stochastic processes in physics and 
chemistry (Elsevier science publishers, 2007). 
E. Arevalo, E.G. Mertens, Y. Gaididei and A.R. Bishop, 
Phys. Rev. E 67, 016610 (2003). 

D.T. Gillespie, Markov processes: an introduction for 
physical scientists (Academic press, 1992). 
D.S. Lemons, An introduction to stochastic processes in 
physics (The John Hopkins university press, 2002). 
M.P. Allen and D.J. Tildsey, Computer simulation of liq- 
uids (Oxford, New York, 1987). 
We tested the range -EswF ~ 10^ — lO'* numerically. 



